New Limits on Interactions between Weakly Interacting Massive Particles and 
Nucleons Obtained with CsI(Tl) Crystal Detectors 
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New limits are presented on the cross section for Weakly Interacting Massive Particle (WIMP) 
nucleon scattering in the KIMS CsI(T^) detector array at the Yangyang Underground Laboratory. 
The exposure used for these results is 24524.3 kg-days. Nuclei recoiling from WIMP interactions are 
identified by a pulse shape discrimination method. A low energy background due to alpha emitters 
on the crystal surfaces is identified and taken into account in the analysis. The detected numbers 
of nuclear recoils are consistent with zero and 90% confidence level upper limits on the WIMP 
interaction rates are set for electron equivalent energies from 3 keV to 11 keV. The 90% upper limit 
of NR event rate for 3.6-5.8 keV corresponding to 2-4 keV in Nal(Ti) is 0.0098 counts/kg/keV/day 
which is below the annual modulation amplitude reported by DAMA. This is incompatible with 
interpretations that enhance the modulation amplitude such as inelastic dark matter models. We 
establish the most stringent cross section limits on spin-dependent WIMP-proton elastic scattering 
for the WIMP masses greater than 20 GeV/c 2 . 

PACS numbers: 95.35.+2, 14.80.Ly 



Astronomical observations have led to the conclusion 
that the majority of the matter in our universe is invis- 
ible, exotic and non-relativistic dark matter How- 
ever, the identity of the dark matter is still not known. 
One possible source are WIMPs, which are candidates 
for particle dark matter that naturally occur in theories 
that extend the standard model of the particle physics 
for reasons independent of the dark matter problem 
There have been a number of experiments that search 
for WIMPs in our galaxy by looking for nuclei recoil- 
ing from WIMP-nucleus scattering 0, U 0, 0^3 ■ To 
date, there are several experiments that interpret their 
results as being possibly due to WIMP signals includ- 
ing DAMA, CoGeNT and CRESST [4-6]. In particular, 
the DAMA results have attracted considerable attention 
since it has reported observations of an annual modu- 
lation of WIMP-like signals with a claimed significance 
of 9 a 4]. This has spurred a continuing debate con- 
cerning the observation of WIMPs that has lasted for 
over a decade. WIMP-nucleon cross sections inferred 
from the DAMA modulation are in conflict with limits 
from other exp eriments that directly measure nuclear re- 
coils @, [9Ml7|- In attempts to reconcile these results, 
various schemes have been suggested, including the in- 
elastic dark matter (iDM) model [19j, in which an ex- 
cited state of the dark matter particle is hypothesized 
and the dominant WIMP-nucleon scattering process in- 
volves a transition to this excited state. Recently, strong 
constraints on the allowed iDM model parameter space 



have been reported H HJ IE | • 

The KIMS (Korea Invisible Matter Search) collabo- 
ration is performing direct searches for WIMPs using 
CsI(T^) detectors in the Yangyang Underground Lab- 
oratory (Y2L). CsI(T^) is a commonly used scintillat- 
ing crystal with 133 Cs and 127 I target elements that 
are sensitive to both spin- independent (SI) and spin- 
dependent (SD) interactions 0]. Furthermore, a pulse 
shape discrimination (PSD) technique makes it possible 
to distinguish nuclear recoil (NR)-induced signals from 
electron recoil (ER)-induced signals on a statistical ba- 
sis. Results from the KIMS experiment based on a four- 
crystal array are reported in Ref. Q- The detector has 
been upgraded to a 3x4 CsI(T^) crystal array with to- 
tal mass 103.4 kg. Each detection module consists of 
a low-background CsI(Ti) crystal (8 x 8 x 30 cm 3 ) Q 
with photomultiplier tubes (PMTs) mounted at each end. 
The data were collected between September 2009 and Au- 
gust 2010. The crystal array is surrounded by a shield 
consisting of: 10 cm of copper; 5 cm of polyethylene; 
15 cm of lead; and 30 cm of liquid-scintillator-loaded 
mineral oil to stop external neutrons and gammas and 
veto cosmic-ray muons. The trigger condition is two or 
more photoelectrons (PEs) in each PMT within a 2/xs 
time window. Amplified PMT signals are digitized by 
400 MHz flash analog-to-digital converters (FADC); the 
total recorded time window for an event is 40 fis, of which 
25 /is is analyzed. Since high energy muon interactions 
in the CsI(T£) crystal can produce a long tail that may 
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last for several tens of milliseconds, we veto events that 
occur less than 50 milliseconds after a muon coincidence 
with the CsI(T£) detector in the off-line analysis. The 
muon coincidence rate is 6-7 per hour and the dead time 
from this veto is negligible. In the energy region below 
10 keV, PMT noise produces a serious background. To 
characterize PMT noise-induced events, a PMT Dummy 
Detector (PDD) module that has the same structure as 
a CsI(T^) detector module with the crystal replaced by 
a clean, transparent and empty acrylic box is included 
in the shielded volume and operated simultaneously with 
the CsI(T^) detector array. The event selection efficiency 
is detector-module dependent and ranges from 20~40% 
at 3 keV to 40-70% at 10 keV. Events that trigger two 
or more detector modules in the array are rejected off- 
line since a WIMP interacts only with the one nuclei in 
a detector: a NR event that is confined to a single detec- 
tor module is the experimental signature for a WIMP- 
nucleus scattering. 

We determine the NR event rate from a PSD analysis. 
To characterize the PSD, we use the quantity we call 
LMT10, which is the logarithm of the mean time of an 
event calculated over a 10 /is interval that starts with the 
first detected PE. Specifically, MT10 (mean time over a 
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where A is the area of the 
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cluster, which is usually equivalent to a single PE, of an 
event. 

We determine the expected LMT10 distribution for ER 
events by irradiating small CsI(T^) crystals (3 x 3 x 1.4 
cm 3 , test crystals 1 & 2) with gamma rays from a 
137 Cs radioactive source and also from Compton scat- 
tering events detected in the array. The test crystals 
are small pieces that have been cut from the same in- 
gots from which crystals used in WIMP search were also 
cut. The expected distributions for NR events were de- 
termined by exposing test crystal 2 to neutrons from a 
Am-Be source. With the higher sensitivity of the current 
exposure, we detected a previously unseen third compo- 
nent with a mean time distribution that is faster than 
that for NR events. Studies show that this originates 
from alpha decays of radioactive isotopes that adhere to 
the crystal surfaces, i.e., surface alph as ( SA), with char- 
actersitics described in detail in Ref. [23(. To determine 
the LMT10 distribution for SA events, we contaminated 
test crystal 1 with Rn progenies, and studied its response 
for events tagged as outgoing alphas from the crystal sur- 
face 23[ . Figure [1] shows a comparison of the three dif- 
ferent reference LMT10 distributions. 

The test-crystal measurements were done at temper- 
atures of (25.4 ± 0.3)°C for SA and (25.3 ± 0.7)°C for 
NR. Throughout the WIMP search data-taking period, 
the crystals were kept in the 20 ~ 21.6°C temperature 
range, with variations that depend on the detector po- 
sition; each detector had rms temperature fluctuations 
of - 0.2°C. The mean values of the ER LMT10 dis- 
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FIG. 1: LMT10 distributions at 3 keV for (a) SA with test 
crystal l,(b) gammas with test crystal 1, (c) neutrons with 
test crystal 2, (d) gammas with test crystal 2 (e) Compton 
scattering events in detector used in the WIMP search. This 
data sample is obtained at (25.2±0.2)°C for comparisons with 
other reference data. 

tributions of the 12 detectors determined in situ have 
an average value of 0.62 with an rms spread of 0.035 for 
59.54 keV gammas from an Am source. Since the tem- 
perature conditions of the reference samples used to de- 
termine the NR and SA distributions are different, some 
adjustment is necessary. Previous studies 24j, [25j have 



shown that the ratio R T = T n jr e is independent of tem- 
perature, where r n is the mean time of NR and r e , the 
mean time of ER. By comparing the LMT10 distribution 
from gamma irradiation on a test crystal and that of the 
crystals used in the WIMP search, we adjust the LMT10 
distribution of NR events and SA events. We assume 
that there are only three components -NR, ER, and SA- 
in the WIMP search data and form the logarithm of the 
likelihood function for each energy bin as: 

Fi = -log(£) 

= - / ] log(fN,jPN,i( x k) + fs,iPsA x k) 
+ 0--fN,i-fs,i)PE,i(Xk)), 



k=l 



where Xk is the LMT10 of the fc th event, the index i 
denotes the i th energy bin, Hi is the number of events 
in the i th energy bin, P/v.j, Pe,i and Ps.i are the prob- 
ability density functions (PDFs) for NR, ER and SA, 
respectively, and /jy,i and fs,i are the NR and SA event 
fractions in the i th energy bin, respectively. 

We determine /jv,i and fs.j using the Bayesian Analy- 
sis Tool (BAT) program [26j with prior PDFs for /jy,, and 
fs.i that are flat between and 1. The Bayesian anal- 
ysis produces posterior PDFs that reflect the degree of 
preference for the parameter values based on the exper- 
imental data. As an example, Fig. [2ja) shows contours 
of the two-dimensional posterior PDFs for /jy,i and fs,i 
for the 6 keV energy bin in one of the detector modules. 
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FIG. 2: a) Two-dimensional plot of /jv,i (horizontal) versus 
fs,i (vertical), b) Projections of fs,i and c) /jv,i for the 6 keV 
bin of detector 9. d) The fitted LMT10 distribution for the 
6 keV bin of detector 9. 

The posterior PDF of each parameter is the projection 
of the two-dimensional PDF onto that parameter's axis 
as shown in Figs.EJb) and (c). To test the sensitivity to 
the choice of priors, we repeated this analysis with the 
flat prior PDF replaced with a Jeffrey's prior [27] . The 
Jeffrey's prior results agree with those with the flat prior 
within 1^2 percent. 

The NR event rate for each detector, determined from 
fjq,i and the event selection efficiency, is shown in Fig. |3] 
The black horizontal bars indicate the 90% confidence 
level (C.L.) upper limits and the red vertical lines denote 
the 68% C.L. intervals. The red horizontal lines mark 
the most probable NR rate values. 
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FIG. 3: Nuclear recoil event rates for all detectors. The black 
horizontal bar indicate 90% C.L. upper limits, the red vertical 
lines denote the 68% C.L. interval, and the red horizontal bars 
the most probable values. 
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FIG. 4: Total nuclear recoil event rates from the combined 
results from nine detectors (without detector 0, 8 and 11). 

The rates of SA events in the 3-11 keV energy range 
averaged over the detector modules are shown in Table HI 
The SA background levels of detectors 0, 8 and 11 are 
about three times as high as the average of the remain- 
ing detectors. For this reason, these three detectors are 
excluded from the average NR event rate determination. 
This reduces the total exposure used in the final analysis 
to 24524.3 kg-days. Figure H shows the 68% C.L. inter- 
vals and 90% C.L. upper limits on the NR event rate 
from the combined PDF of the remaining nine detector 
modules. These limits include systematic effects from un- 
certainties in the LMT10 PDFs for each event type and 
their crystal-to-crystal deviations. No significant excess 
of NR events is observed. 

Assuming the standard halo model fUJ, we translate 
these measurements to the cross section limits for WIMP- 
nucleon SI interactions and WIMP-proton SD interac- 
tions that are presented in the left and right panels of 
Fig. 03 respectively. The limits shown in the two figures 
are about one order-of-magnitude more stringent than 
the previous KIMS results [7] , due primarily to the iden- 
tification of the SA background component and the larger 
exposure. The WIMP-proton SD interaction cross sec- 
tion limits are the most stringent to date. 

The NR event rate limits also have important im- 
plications for the interpretation of the DAMA annual 
modulation signal, which has an amplitude of 0.0183 ± 
0.0022 counts/day/kg/keV in the 2-4 keV energy range 
in NaI(T£) scintillators Considering the different 

quenching factors of Iodine for Nal(T^) and CsI(T£) [3, 
[29j, the 2-4 keV DAMA energy range corresponds to 
3.6-5.8 keV in KIMS, which is included in the first 
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0.203 ± 0.026 
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0.087 ± 0.021 
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0.071 ± 0.017 
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0.076 ± 0.025 
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0.066 ± 0.020 
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0.238 ± 0.025 


3 


0.089 ± 0.024 
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0.123 ± 0.025 
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0.039 ± 0.020 


10 


0.014 ± 0.026 
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0.064 ± 0.018 


11 


0.205 ± 0.024 



TABLE I: SA contamination level for each detector averaged 
over the 3-10 keV energy bins. 
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FIG. 5: The 90 % exclusion limits on (Left) SI WIMP-nucleon and (Right) SD WIMP-proton cross sections. In both plots 
DAMA results interpreted by Savage et al . 11711 are used (3a contours are drawn). The SI plot includes NAIAD [9|], CRESST- 
II [3, EDELWEISS-II [n]], ZEPLIN-III [lj XENON100 QJ] and CDMS limits. The SD plot includes PICASSO HH 
and COUPP [3 limits. 



three bins in Fig. 2] Our 90% C.L. upper limit on 
the NR event rate in the 3.6-5.8 keV energy range 
is 0.0098 counts/day/kg/keV, which is well below the 
DAMA signal amplitude. Therefore, any scenario involv- 
ing Iodine as the target, such as the iDM model, is incom- 
patible with our limits. As an example, the parameter 
space allowed for DAMA in the iDM model and our ex- 
clusion limits for a WIMP of mass 70 GeV are presented 
in Fig. [5] An alternative iDM interpretation considers 
Thallium, which is present at the 10 -3 level in both the 
DAMA and KIMs detectors 22], as the dominant target, 
can be addressed by our results. We estimate the quench- 
ing factors for Thallium in NaI(T£) and CsI(T£) using 

a semi-empirical calculation 30] and find ^p 1 l^fS 

where Qcili NaI * s tne quenching factors of CsI(T^) and 
Nal(T^) for iodine and thallium ions. This indicates that 
the corresponding energy range in KIMS for Thallium is 
about the same as that for Iodine. Therefore our conclu- 
sion does not change when Thallium is considered as the 
dominant target. 

In conclusion, we report improved limits for WIMP- 
nucleon cross sections using a data sample collected with 
a 103.4 kg CsI(T£) scintillator detector array with a total 
exposure of 24524.3 kg-days. We identified and charac- 
terized a low energy background due to a contamination 
of alpha emitters on the surfaces of the crystals and incor- 
porated it into the PSD analysis. No significant signals 
for NR events are observed and we determine 90% C.L. 
upper limits on NR event rates, and improved limits on 
WIMP-nucleon cross sections, including the most strin- 
gent limits to date on WIMP-proton SD scattering. The 
NR event rate upper limit is below the DAM A/LIBRA 
annual modulation amplitude in the corresponding en- 
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FIG. 6: The allowed parameter space for DAMA/LIBRA [l<J 
and the limits reported here for a 70 GeV WIMP mass in 
iDM model. 5 is the mass split between the ground and ex- 
cited states of the WIMP. The astronomical parameters from 
Ref. [l9| ] are used. 



ergy region, disfavoring iDM model interpretations. 
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